Nonlinear stress and fluctuation dynamics of sheared disordered wet foam 
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Sheared wet foam, which stores elastic energy in bubble deformations, relaxes stress through 
bubble rearrangements. The intermittency of bubble rearrangements in foam leads to effectively 
stochastic drops in stress that are followed by periods of elastic increase. We investigate global 
characteristics of highly disordered foams over three decades of strain rate and almost two decades 
of system size. We characterize the behavior using a range of measures: average stress, distribution 
of stress drops, rate of stress drops, and a normalized fluctuation intensity. There is essentially no 
dependence on system size. As a function of strain rate, there is a change in behavior around shear 
rates of 0.07 s~\ 

PACS numbers: 83.80.Iz,82.70.-y 



I. INTRODUCTION 

One of the potentially exciting features of driven, com- 
plex fluids is the possible existence of an "effective" tem- 
perature [ij, |2| • Examples of systems for which an effec- 
tive temperature may prove a useful idea include foams, 
emulsion, granular materials, and colloidal glasses (for 
example, see Ref. |3j, and references therein). Theoret- 
ical studies of effective temperatures using the bubble 
model of foams pj and the "standard model" for super 
cooled liquid (a binary Lennard- Jones mixture) 0] pro- 
vide strong motivation for experimental studies of effec- 
tive temperature. Understanding the nature of fluctua- 
tions in these systems is a key step toward developing 
an understanding of any concept of effective tempera- 
ture. In this paper, we focus on fluctuations in a sheared, 
two-dimensional foam system: bubble rafts 0, Q . Un- 
der shear, initially jammed foam exhibits elastic behav- 
ior until the yield strain, followed by intermittent "stress 
drops" as the bubbles undergo nonlinear topological re- 
arrangements (for a review of foams see 0, @| ) ■ The 
fluctuations in stress, and other quantities in the system, 
are reminiscent of thermal fluctuations and motivate defi- 
nitions of effective temperature. Previous work with bub- 
ble rafts characterized the statistical distribution of stress 
drops for a single system size and a small range of strain 
rates [5j. The results were in excellent agreement with 
simulations of the bubble model P. Hollri| . In this paper, 
we will report on results for a greater range of system size 
and strain rate. Additionally, we will report on a number 
of measures other than the distribution of stress drops. 
In addition to providing a starting point for studies of 
effective temperatures, this work provides detailed tests 
of competing models of flowing, two-dimensional foams. 

There are a number of different models of flowing foam. 
They all make qualitatively similar predictions regarding 
the behavior of the stress as a function of strain rate. 
For small rates-of-strain, there is an initial elastic region. 



At a critical value of the stress (yield stress) or strain 
(yield strain), the foam begins to flow. The flow in this 
region is intermittent, with periods of increase followed 
by sudden, irregular stress releases, referred to as stress 
drops. Where the models differ is the details of the dis- 
tribution of stress drops. These differences are the result 
of different assumptions concerning the source of dissipa- 
tion and dryness of the foam. As this work focuses on 
two-dimensional foams, the "dryness" is characterized by 
the area fraction of gas, 0. For 0=1, foam is perfectly 
dry, and the bubbles are all polygons. For <p < 0.84, the 
foam "melts" into a froth of exclusively circular bubbles. 
Foams near this limit are referred to as "wet" . 

Four main models of two-dimensional foam are: the 
vertex model 0,0,0; the quasi-static model 0,0]; 
the extended q-Potts model [l7|; and the bubble model 
0, ^3 ■ The vertex model is not particularly relevant 
to our system, as it models dry foam. The quasi-static 
model is special because it docs not contain any dissi- 
pation. It deals with wet foams. The results from this 
model suggest a power-law distribution for the proba- 
bility of stress drops of a certain size occurring (lEL Il6| . 
This result is based mainly on measuring the distribution 
of Tl events. Tl events are nonlinear neighbor switch- 
ing events between four bubbles. The extended q-Potts 
model includes dissipation without making any specific 
assumptions about the dissipation. The work in Ref. 0] 
focuses on the dry foam limit, but the model can treat 
wet foams. These simulations suggest that a quasi-static 
limit does not exist for foam. In other words, as the 
strain rate is continually decreased, the properties of the 
flow continue to change. Also, they report a transition, 
as a function of disorder, from a viscoelastic solid to a 
viscoelastic fluid. The signature of this change was the 
fact that sufficiently disordered foams displayed no yield 
strain. Instead, such foams flowed immediately upon the 
application of shear. Finally, this work suggests that 
the distribution of energy drops in sheared foam obeys a 
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power-law. However, the distribution of Tl events does 
not. The bubble model is applicable to wet foams, and 
its main prediction with regard to stress drops is a power- 
law distribution for small stress drops with an exponen- 
tial cutoff at larger stress drop magnitudes [3, Hi! llil ■ 
In contrast to the q-Potts model, the bubble model does 
predict a quasi-static limit for flowing foams. 

Experiments with foam have yielded various results. 
Indirect studies of three-dimensional foam using diffu- 
sive wave spectroscopy suggest a cutoff to the non- 
linear rearrangements jl8j . Studies of Tl events us- 
ing two-dimensional foams agree with these results |l9j . 
These experiments used the gas-liquid coexistence region 
of Langmuir monolayers to make truly two-dimensional 
foams. Experiments using quasi-static shear of a single 
layer of bubbles between glass plates suggest that there 
may be system-wide events, suggesting a possible power- 
law behavior [2(j ■ Work with bubble rafts (a single layer 
of bubbles on the surface of water) showed that the distri- 
bution of stress drops for a system of approximately 900 
bubbles exhibited the exponential cutoff predicted by the 
bubble model |5| . In order to further test the agreement 
with the bubble model, we have extended the work in 
Ref. Q to a wide range of system sizes and strain rates. 
Furthermore, in addition to reporting on the distribu- 
tion of stress drops, we also measure the rate of stress 
drops, the average stress (which gives the viscosity), and 
the normalized stress fluctuation. We establish that the 
cutoff in stress drops is not a finite size effect. Also, we 
show that for this system there is a well-defined quasi- 
static limit, as predicted by the bubble model. Finally, 
we will discuss the apparent dependence of some of the 
properties on system size and strain rate. The rest of the 
paper is organized as follows. Section II describes the ex- 
perimental techniques. Section III presents the results. 
Section IV discusses the results in context of the various 
models of foams. 



II. EXPERIMENTAL DETAILS 

A standard bubble raft Q provides a nice model sys- 
tem for studying two-dimensional foams. Bubble rafts 
consist of a layer of bubbles floating on the surface of 
water. The motion of the bubbles is essentially all in the 
plane defined by the water surface. However, it should 
be noted that the system is not an ideal, two-dimensional 
system, as some motion is possible perpendicular to the 
surface. However, in all of our experiments, the bub- 
bles were monitored with video cameras, and no motion 
was observable perpendicular to the water surface. The 
bubble raft was produced by flowing regulated nitrogen 
gas through a hypodermic needle into a homogeneous 
solution of 82.0% by volume glycerine, 14.5% by volume 
deionized water, 1.50% by volume tricthanolamine, and 
2.00% by volume oleic acid. The bubble size was depen- 
dent on the nitrogen flow rate, which we varied using a 
needle valve. The bubble diameter ranged from 2 to 6 



mm, with most bubbles in the 3 to 4 mm range. The re- 
sulting bubbles were spooned into a cylindrical Couette 
viscometer described in detail in Ref. |2lJ. This pro- 
duced a two-dimensional, wet foam on a homogeneous 
liquid substrate of 80% by volume deionized water, 15% 
by volume glycerine, and 5.0% by volume Miracle Bub- 
bles (Imperial Toy Corp.) The Couette viscometer con- 
sists of a shallow dish that contains the liquid substrate. 
Two concentric Teflon rings are placed vertically in the 
dish. The outer ring consists of twelve segmented pieces 
and has an adjustable radius. The inner ring, or rotor, 
has a radius r = 4.0 cm and was suspended by a wire to 
form a torsion pendulum. Polypropylene balls with a 4 
mm diameter were epoxied to the circumference of the 
inner rotor to prevent the innermost row of bubbles from 
slipping. The outermost row of bubbles was monitored, 
and no slip of bubbles on the outer barrier was observed. 
Figure ^ is a top view of a typical bubble raft in our 
apparatus. Portions of the two Teflon rings are visible. 

To shear the foam, the outer Teflon barrier was ro- 
tated at a constant angular velocity in the range 0.0005 
rad/s to 0.5 rad/s. The torque, r, on the inner rotor 
was monitored by recording the angular position of the 
inner rotor (and thus the angular displacement of the 
torsion wire) twice per second. The angular position was 
measured using magnetic flux techniques described in de- 
tail in Ref. |2l|. This data was recorded by a 12-bit 
analog-to-digital converter and stored on a PC. The tan- 
gential stress, a, on the inner rotor due to the foam is 
given by a = r/(27rr 2 ). During periods without any re- 
arrangements, the fluctuations in the stress were at the 
level of one bit, corresponding to changes in stress of 
2 x 10~ 3 dyne/cm. Therefore, when computing stress 
drops, we filtered the data to eliminate any changes in 
stress of ±2 x 10 -3 dyne/cm. 

The bubble raft was constructed by placing the ap- 
proximate number of desired bubbles in the trough with 
the outer barrier set to a large radius. It is important 
to note that the bubbles exhibited a strong attraction 
to each other. This is a microscopic detail that is not 
included in any of the models discussed in the introduc- 
tion. Typically, two-dimensional foams are characterized 
by their gas-area fraction, which is the ratio of the area 
filled by gas to the total area. Because the bubbles actu- 
ally exist in three dimensions for a bubble raft, the fluid 
walls (and the cross-sectional area of a bubble) are height 
dependent. This complicates the definition of gas-area 
fraction. Therefore, we used a functional definition for 
gas-area fraction based on the images of the bubble raft. 
We defined the area of gas to be the black regions within 
bubbles in an image and maintained constant lighting 
conditions so that this definition was consistent from run 
to run. The outer barrier was compressed until the de- 
sired bubble density was achieved. It should be noted 
that this resulted in a variation in the initial shear stress 
of the bubble raft that did not relax significantly on the 
time-scale of the experiments. Therefore, due to the fi- 
nite lifetime of the raft, the experiments were carried 
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out with this initial pre-stress present. Both the total 
number of bubbles and average gas area fraction were 
determined from images of a large section of the trough, 
assuming an essentially uniform distribution of bubbles 
throughout the trough. For all of the data reported here, 
the gas area fraction was approximately 0.95. 

The stability of the bubbles was enhanced by cool- 
ing the fluid substrate to 5 °C. Also, a glass cover was 
placed over the bubbles. The cover helped reduce evap- 
oration and was not in contact with the bubbles. The 
entire apparatus was contained in a cabinet. The cabi- 
net reduced the air flow around the apparatus, and a hu- 
midifier placed within the cabinet helped to extend the 
lifetime of the bubbles. The bubbles in the bubble raft 
did not exhibit any substantial coarsening with time. In- 
stead, the raft tended to suffer catastrophic failure after 
approximately two hours due to a significant number of 
bubbles popping. Presumably, this was due primarily to 
the loss of fluid in the bubble walls from drainage into the 
fluid substrate and/or evaporation. Because the bubble 
raft did not coarsen significantly, there is no competition 
between coarsening and shear induced rearrangements. 



III. EXPERIMENTAL RESULTS 

Figure El shows the typical behavior of the stress ver- 
sus strain for strains above the yield strain. The irregu- 
lar behavior of the stress during flow is apparent. This 
behavior can be characterized by considering the distri- 
bution of stress drops. This is shown in Fig. [3] for a 
system with 1550 bubbles for various strain rates (7). 
All results are given in terms of normalized stress drops 
(Act = 5<J / Umax), where the stress drops 5a for each run 
are normalized by the maximum stress {a max ) for that 
run. A bin size for Act of 2.5 x 10~ 3 was used in plotting 
the distribution, with the probability of a stress drop of 
size Act (P(Act)) defined as the number of drops within 
each bin divided by the total number of drops in the run. 
The solid line is a guide to the eye and has a slope of 
-0.8. The distribution is consistent with a power law for 
small stress drops with an exponential cutoff. Because of 
the cutoff, there is a a well-defined average stress drop, 
< Ac >. The P(Act) for different system sizes are qual- 
itatively the same, a power-law for small Act, with an 
exponential cutoff. In order to look for quantitative dif- 
ferences as a function of system size, we considered the 
behavior of < Act >. 

The average stress drop is shown as a function of strain 
rate for different system sizes in Fig. E| The smallest 
system consisted of 1.6 x 10 3 bubbles and is given by 
the squares. For this system, there is essentially no de- 
pendence of < Act > on strain rate. For the systems 
with more than 5.6 x 10 3 bubbles (all of the other sys- 
tems that we studied), there does appear to be a weak 
strain rate dependence. One observes an increase in the 
average stress drop with strain rate, until 7 « 0.07 s . 
Above this value, the average stress-drop is independent 



of strain rate. The implications of this will be discussion 
in Sec. IV. 

One important result is that there is no increase in 
< Act > with system size. This effectively rules out 
system size as the source of the cutoff. At the lowest 
strain rates, < Act > is lower for the larger systems. A 
possible reason for this behavior is discussed in Sec. IV. 

The relation between average stress and strain rate is 
shown in Fig. |SJ The maximum stress displays a similar 
dependence. Note the knee of the curve, at 7 « 0.07 s _1 . 
A line is drawn as a guide to the eye, with a slope of 1/3. 
The overall curve is consistent with a Herschel-Bulkley 
model of viscosity where the stress is given by ct = A + 
B7™ [23. Here 7 is the strain rate and A and B are 
constants. The variation in average stress in Fig. [5] is 
most likely due to the variation in the initial stress from 
run to run discussed in Sec. II. Presumably, if sufficient 
aging of the system were possible before each run, this 
variation would be reduced. 

An alternate way to view the same data is to con- 
sider 77 = < ct > /j versus 7, shown in Fig. [B] This 
is the steady-state viscosity, taking care in computing 7 
at the inner cylinder [2^J . The solid line with slope -1 
and dashed line with slope -2/3 are guides to the eye 
and clearly illustrate the Herschel-Bulkley behavior of 
the bubble raft, with an exponent n = 1/3. This behav- 
ior is consistent with the shear-thinning velocity profile 
reported in Ref. |5] for 7 = 0.062 s . An open question 
is the behavior of the velocity profile at extremely low 
strain rates, where the average stress is essentially inde- 
pendent of strain rate. One definitely observes bubble 
motions throughout the bulk of the system. However, 
at these low shear rates, most of the time the system is 
undergoing a linear increase of the stress, and only occa- 
sionally is there a stress drop. Initial measurements of the 
flow during such an increase in the stress are consistent 
with a linear profile of the velocity. However, detailed 
measurements in this regime will be conducted in the fu- 
ture to determine what the long-time average (one that 
includes many stress drops) of the flow profile is. This 
is important given measurements of velocity profiles at 
low shear rates that report an exponential decay of the 
velocity for a similar system |23|. For our system, one 
issue is whether or not the water substrate dragged the 
bubbles. We made a number of measurements where the 
outer cylinder was rotated, but the bubble raft was not in 
contact with the outer cylinder. This was accomplished 
by removing the approximately the outer three rows of 
bubbles. Under these conditions, no flow of the bubble 
raft was observed, and no measurable stress was trans- 
mitted to the inner rotor. This provides strong evidence 
that the underlying water does not "drag" the bubble 
raft. 

The dependence of the maximum stress on strain rate 
suggests the existence of a quasi-static limit. Below a 
strain rate of approximately 0.07 s , the stress is essen- 
tially independent of the rate of strain. This was checked 
by considering the number of stress drops per unit strain 
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(S). This is plotted in Fig. As with the maximum 
stress, S approaches a constant below values of the strain 
rate of approximately 0.07 s _1 . 

In addition to the distribution of stress drops, we also 
characterized the intensity of fluctuations around the 
mean. We defined the fluctuation intensity, T, as the 
standard deviation of the stress (after the yield stress) for 
a given run, expressed as a fraction of the mean stress: 



r = 



\ 



1 N 

-Y 

i=l 



< a > 



< a > 



(1) 



where the sum is over the measured values of stress and 
N is the number of data points for a given experimental 
run. Given the existence of the pre-stress, the normaliza- 
tion by the average stress allows for a better comparison 
between different systems. Figure [H] shows the results 
for the fluctuation intensity (r) as a function of 7 for 
the different system sizes. There is significant scatter in 
the data; however, there is a clear trend of decreasing 
r as 7 increases. One consequence of this is a correla- 
tion between T and the < ct >. This is illustrated in 
Fig. El where the data for the different system sizes are 
combined into a single plot. Again, there is significant 
scatter in the data, but the trend is obvious. 



IV. SUMMARY 

The results presented here provide strong evidence that 
the bubble model provides an accurate description of the 
shear behavior of a bubble raft. To the extent that bubble 
rafts are equivalent to foam, the bubble model would also 
describe a two-dimensional flowing foam. The stress drop 
distribution, the average stress as a function of strain 
rate (the steady state viscosity), and the rate of stress 
drops, S, are all consistent with the bubble model [Tlj . 
This agreement is despite the fact that the bubble rafts 
studied here are strongly attractive, a feature that is not 
explicitly in the model. One interesting result is that 
both the bubble model and the bubble rafts are well- 
described as a Herschel-Bulkley fluid with an exponent 
of 1/3 24]. It remains to be seen if this exponent is a 
generic feature of the model, and of similar models, or 
if there is something specific to the parameters used in 
Ref. [2^|. For example, the exponent may depend on 
various characteristics of the foam, such as the gas-area 
fraction. This dependence also needs to be tested for the 
bubble rafts. 

Comparison of our studies with both the quasi-static 
model 0,^1 an d the quasi-static experiments [2(| raises 
an interesting question: is there a fundamental difference 
between slow but steady shear, and true quasi-static mo- 
tions? The disagreement between our results and the 
quasi-static experiments suggests that such a difference 
may exist. However, it is also possible that the discrep- 
ancies are due to comparing direct measurements of the 
stress drops with sizes of spatial rearrangements. Future 



work with our system will look at both the issue of quasi- 
static steps versus steady shear and the spatial extent of 
rearrangements. 

With regard to the extended q-Potts model, this work 
raises some important questions. Two clear predictions 
of that model are: (1) there is no quasi-static limit; and 
(2) a sufficiently disordered foam no longer has a yield 
strain |17| . Neither behavior was observed in our experi- 
ments. At this point, one would need to do further work 
to determine if there was something fundamentally miss- 
ing from the q-Potts model that results in this disagree- 
ment. The other possibility is that our foams were either 
not sufficiently disordered to be accurately described by 
the q-Potts model or they were not sufficiently dry, as 
the simulations in Ref. Qjfl were for dry foam. There- 
fore, future experiments will focus on the role of disorder 
and the wetness of the foam. 

Though not conclusive, the behavior of the average 
stress drop as a function of system size and strain rate, 
shown in Fig. 01 suggests some interesting behavior. The 
large decrease in < Act > as a function of system size is 
surprising. One possible explanation involves the spatial 
correlations between bubble rearrangements that pro- 
duce the stress drops. As the system size increases, there 
are more spatial locations at which rearrangements can 
occur. For low enough strain rates, there will be an in- 
termediate range of systems size for which this increases 
the probability of isolated small stress drops occurring. 
Once one region slips, enough stress is relieved that the 
other regions do not rearrange until a sufficiently later 
time that they are recorded as a new stress drop. Such 
dynamics would result in a decrease in the average stress 
drop with system size. Eventually, as the system size 
increases even more, this behavior should "smooth" out 
the dynamics, as small stress drops occur almost con- 
tinuously. Presumably, this happens in large, three- 
dimensional samples. On the other hand, for these in- 
termediate size systems, as the strain rate is increased, 
the stress releases occur closer together. This increases 
the likelihood of multiple small events in different spatial 
locations combining to form larger stress drops. There- 
fore, one observes an increase in the average stress drop 
as a function of strain rate. Eventually, as the system 
crosses over to more fluid-like behavior, there is again a 
"smoothing" of the dynamics. In this regime, the average 
stress drop becomes independent of strain rate. Clearly, 
more work is needed, both in the experiments and simu- 
lations, to test these ideas. In particular, they highlight 
the importance of measuring both spatial correlations be- 
tween rearrangement events and the correlation between 
the rearrangements and the stress drops. 

Finally, it interesting that the crossover to "smoother," 
more fluid-like behavior, as a function of strain rate, is 
evident in both the measurement of T and < Act > . How- 
ever, the two measures reveal slightly different behavior. 
As one increases the strain rate, T decreases monotoni- 
cally (see Fig. [HJ- However, < Act > becomes indepen- 
dent of 7 at the higher strain rates (see Fig. The 
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two results are not inconsistent, as F measures the size 
of fluctuations from the mean, while < Act > measures 
the average size of the changes in stress that produce 
these fluctuations. Again, understanding the spatial dis- 
tribution of bubble rearrangements will probably be an 
important step in fully understanding the dependence of 
these two measures of the fluctuations on strain rate. 
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FIG. 1: Image of a section of a typical bubble raft. A portion 
of the inner rotor is visible in the lower left corner. A portion 
of the outer barrier is visible in the upper right corner. This 
particular raft had approximately 10 4 bubbles. The scale bar 
is 1 cm. 
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FIG. 2: Typical stress response of a sheared wet foam il- 
lustrating the intermittent stress drops. The rate of strain is 
0.014 s~\ and the number of bubbles is 1.6 x 10 3 . 
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FIG. 3: Probability distribution for stress drops (P(A<t)) as a 
function of the magnitude of the stress drop (Act = Sa/tJmax)- 
For each separate run, the stress drops (Sa) have been nor- 
malized by the maximum stress {o max ). The symbols are 
different strain rates: 2.7 x 10~ 3 s" 1 (v); 1-4 x 10~ 3 s^ 1 (A); 
7 x 10~ 3 s^ 1 (<); 1.4 x 10~ 2 s" 1 (o); 1.3 x 10" 1 s" 1 (0); 
2.7 x 10 _1 s _1 (□). The solid line is a guide to the eye and 
has 
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FIG. 4: Plot of the average size of the stress drop (< Act >) as 
a function of the strain rate for different system sizes: 1.6 x 10 3 
(□); 5.6 x 10 3 (o); 9.2 x 10 3 (A); 1.5 x 10 4 (v); 2.0 x 10 4 (0); 
2.6 x 10 4 (<). The solid vertical line is at j = 0.066 s" 1 . 
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FIG. 5: Plot of the average stress versus the strain rate. The 
dark line is a guide to the eye, and has slope 1/3. 



1 

0.001 0.010 0.100 1.000 

strain rate (s 1 ) 

FIG. 6: Plot of average stress divided by the rate of strain FIG. 9: The fluctuation intensity (F) versus the mean stress. 
(rj) versus the rate of strain. The solid line has a slope of -1 The dark line is a guide to the eye, and has slope -1.3. 
and the dashed line has a slope of -2/3. The two lines cross 
at 7 = 0.066 s" 1 . 
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FIG. 7: The number of stress drops per unit strain (S) as 
a function of strain rate. The plateau at low rates of strain 
suggests a quasi-static limit is reached. 
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FIG. 8: . Fluctuation intensity (F) versus strain rate for a 
range of system sizes: 1.6 x 10 3 (□); 5.6 x 10 3 (o); 9.2 x 10 3 
(A); 1.5 x 10 4 (v); 2.0 x 10 4 (0); 2.6 x 10 4 (<). 



